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Abstract 

We expose a grave drawback hidden in the standard model of nn transitions 
in a medium. Correcting it, one improves a limit on the free-space nfi oscillation 
time T n n by 31 orders of magnitude. 
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Recently, the nn transitions in a medium beyond the potential model have been 
considered [1]. For lower limit r m i n on the free-space nn oscillation time we have got 
T~nn > Tmin = 4.7 • 10 31 y, which increases the previous one by 31 orders of magnitude. 
Certainly, the value 10 31 provokes distrust. In this letter we explain in very simple 
terms the origin of this perceptible (and expected by us) disagreement. 

In the standard approach (labelled bellow as potential model) the nn transitions in 
a medium are described by Schrodinger equations 

(id t + V 2 /2m - U n )n(x) = en(x), (id t + W 2 /2m - U n )n{x) = en(x). (1) 

Here e = l/r nn is a small parameter [2]; U n and U fl are the self-consistent neutron 
potential and n-nucleus optical potential respectively. For U n = const, and U n = const. 
in the lowest order on e the probability of the process is 

W pot (t) = 2ImTii(t), T u {t) = t(e/5U) 2 [l - iSUt - exp(-i5Ut)], (2) 

where 



5U = U n -U n , U n = ReUn-iT/2, (3) 

T ~ lOOMeV is the annihilation width of n-nucleus state. The limit for r nfl is obtained 
from the inequality W pot (T) < 1, where T ~ 10 31 y is the experimental bound for the 
nuclear annihilation lifetime. Priviously, we adduced a physical arguments that the 
model (1) is invalid[l,3]. Now we will pinpoint the origin of error exactly. 

W pot describes the nn-transition, annihilation. Let us consider the second stage - 
annihilation decay of n-medium state. The wave function of initial state obeys equation 

H <f> = e n $ : H = -V 2 /2m + U n . (4) 
In r = the interaction 5U is turned on. We have 

— = (H + SU)y, (5) 
\l/(0) = $. The projection to the initial state is 

< $ | *( T ) >= Uu(t) = exp(-i5Ur), (6) 

where U(t) — 1 + iT n {r) is an evolution operator, T n (r) is a transition matrix. The 
decay probability Wn(r) can be obtained by two ways. 
1. First one gives correct result: 
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W n (r) = 1- | Uu(t) | 2 = 1 - e~ v \ (7) 

2. Taking into account that unitarity condition UU + = 1 must be fulfilled for 
[/-matrix of any process we have 

W n {r) = 2/m7™(r). (8) 
Substituting T-™ = i[l — exp(—i5Ur)] from Eq.(6), one obtains 



Wjr) 



1 - e' Fr/2 cos(rRe5U)] , (9) 



which strongly differs from (7). When r >> 1/T ~ 10 24 s, Eq.(9) gives 

W n = 2. (10) 

The probability nonconservation was to be expected because U n is non-Hermitian and 
unitarity condition is strongly violated: UU + = e~ Tr — > 0. As we will see bellow such is 
indeed the case in the calculation of (2). In the slight absorption region Yr ~| 5Ut \« 
1 the probability nonconservation is small: UU + — 1 — Tt m 1. As a consequence of 
this Eq.(9) gives result identical to Eq.(7) one: W n = Tt. 

Therefore, the matrix element ImT^ constructed with the help of non-Hermitian 
operator U n has a sense only when W n — Tr << 1, or, what is the same 

2M™(r) « I. (11) 

When Ft > 1, expression for ImT£ (i.e.,Eq. (9)) must be replaced by (7), which is 
what we made in Ref.[l]. The alternative is to calculate ImT£ beyond the potential 
model. In any case the equality 2ImT^{r) = 1 — exp(— Ft) is a test for T£ calculation. 

Let us return to Eq.(2). For the whole process probability we know full well that 
Wp t « 1. The probability nonconservation is also small UU + = 1 — 4e 2 t/r ss 1. 
Thus, as is easy to verify, the both methods give the same result (2). However, the 
matrix element T {i itself, and with it the process probability W pot are in error because 
they are calculated by means of the erroneous "decay law" (9) in the region Ft >> 1. 
In the potential approach this fact is strongly hidden. In the approach with finite time 
interval [1] it is obvious. 

Really, instead of (2) we have[l] 

W{t) = W a (t) + W b (t), W a (t) = e 2 t 2 (12) 
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W\t) = -e 2 j dt a I 
Jo Jo 



dtpWfiita-tp), 



(13) 



where W a ,W b are the contributions of Fig. 2a and 2b of Ref.[l] respectively, Wn(r) is 
the probability of n-nucleus decay in time r = t a — tp. Substituting Eq.(9) in (13), one 
obtains 



Consequently, the potential model result is really obtained by means of "decay law" 
(9) and in doing so the main contribution gives the range where W n ~ 2. (The latter 
is apparent from (13), where t — T — 6.5 • 10 31 y[4].) Certainly, this can be understood 
without going beyond the potential model framework. Indeed, solving the system (1) 
by method of Green functions we get 



Now we have to ask himself: What is used as matrix element of annihilation decay Tf { 
in (15)? Eq.(lO) provides the answer. 

The conclusion that potential model is inapplicable to the problem under study[l] 
means as follows. Obviously from the point of view of microscopic theory the potential 
approach is in principal invalid. Certainly any characteristic can be parametrized. 
Eq.(7) is an example of useful parametrization. However, instead of it erroneous Eq.(9) 
is used. This inevitably arises in solving Eqs.(l), because 2ImTa is represented through 
2ImTZ , other than 1- | U u {t) \ 2 . 

Our result is obtained by substituting of exponential decay law (7) in (13): 



Moreover, in Eq.(13) one can put W n = const < 1. (We would like to stress this 
circumstance.) Then W > e 2 t 2 /2. In this case only unitarity condition was employed 
in deriving of W(t). 

It remains to see the reason of enormous quantitative disagreement between the 
our and potential model results. The formal answer is that the overall factor 2 in 
Eq.(9), compared to Eq.(7), leads to the full cancellation of the e 2 t 2 terms in Eq.(14). 
The strong result sensitivity was to be expected. Really, the S- matrix amplitude M s , 
corresponding to nn transition, annihilation diverges (see Fig.l and Eq.(5) of Ref.[l]): 



W(t) = e 2 t 2 - e 2 t 2 + W pot (t) = W pot (t). 



(14) 




(15) 




(16) 
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(17) 



M s 



= e 



P 2 /2m - U, 







n 



where M is the annihilation amplitude. This is infrared singularities conditioned by 
zero momentum transfer in the e- vertex. It is easy to understand that M s ~ 1/0 for 
any bound state wave function of neutron (i.e., for any nuclear model). On the other 
hand from Eqs.(l) it is clear that in the potential model the energy is not conserved 
and becomes complex in the e- vertex Ma — ► Ma + SU (Ma is the nuclear mass). The 
corresponding antineutron Green function is 



SU = is the peculiar point of M s . So M s is extremely sensitive to SU ; the result W(t) 
is extremely sensitive to W n {r) obtained by means of SU. (Usually, the ^[/-dependence 
of G is masked by q: G~ l = (e n — q ) — (p n — q) 2 /2m — U n — SU. We deal with 2-tail 
and q = 0.) 

The nn transition, annihilation (two-step nuclear decay) and particles motion in 
the classical fields are the different problems. Describing the first one by Eqs.(l) we 
understand that this is an effective procedure. From formal standpoint in the first 
and second cases the potentials are complex and real respectively. (Unfortunately, 
sometimes the literal analogy between these problems is drawn[5].) If we try to employ 
the potential description (1) as effective one we must invoke an additional reasoning 
for verification. Otherwise, we can get (10). 

Therefore, the potential approach is inapplicable to the problem under study nei- 
ther from microscopic theory standpoint nor from phenomenological parametrization 
one, because condition (11) is not fulfilled. More expended conclusion is as follows. 
The energy gap SU ^ is responsible for dramatic process suppression. It is inevitably 
exists in the potential model. However, on the reasons presented above, the potential 
approach is in principle invalid. This explains essentially different functional struc- 
ture of the results: W(T)/W pot (T) ~ FT. If it is remembered that T ~ 10 31 y the 
quantitative distinction becomes clear as well. 



G = l/(e n - p^/2m - U n ) = l/SU. 



(18) 
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